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Abstract. We give asymptotic analysis for probability of absorbtion 
P(ro < T) on the interval [0,r], where to — inf{t : Xt — 0} and Xt is a 
nonnegative diffusion process relative to Brownian motion Bt, 

dXt = ^iXtdt + aX^dBt. 
Xo = K >0 

Diffusion parameter ax'' , 7 G [^,1) is not Lipschitz continuous and 
assures P(ro > T) > 0. Our main result: 

where Mr = /o^cr(l - 7)e-<i-T''"'dB^. Moreover we describe the most 
likely path to absorbtion of the normed process ^ for K ^ 00. 

1. Introduction 

In this paper, we analyze the Constant Elasticity of Variance Model 
(CEV), introduced by Cox 1996, [1] and appUed to Option Pricing (see 
e.g. Delbaen, Shirakawa [2] and Lu, Hsu [9]). This model is given by the 
Ito equation with respect to a standard Brownian motion Bt and a positive 
initial condition Xq = K > 0, 

dXt = fiXtdt + aX^dBt, (1.1) 

where /i, o" 7^ are arbitrary constants and 7 € 1). For 7 = ^, this model 
is known as CIR (Cox, Ingersol and Ross) model. The diffusion coefficient 
ax'^ is only Holder continuous, yet the Ito equation (1.1) has unique strong 
solution^. In contrast to Black-Scholes model (7 = 1) with Xt > for 
any t > 0, for CEV model the process Xt is absorbed in zero at the time 
To = inf{t : Xt = 0} with P{tq < 00) > which can be interpreted as time 
of ruin. 

In a proposed asymptotic analysis, as X ^ 00, a crucial role plays the 
normed process xf = being the unique solution of the following Ito 
equation 

dxf = ^ixfdt + -^{xfydBt, (1.2) 

subject to the initial condition Xq = 1, and a small diffusion parameter 
j^f^^ . We emphasize that the process xj^ inherits the ruin time tq. 

The assumption 7 < 1 implies that the diffusion in (1.2) has a small 
diffusion coefficient. This enables us to find a rough lower bound of P(ro < 
T) for any K > (see Remark 2). With K — > 00, this lower bound is best 



"'^Delbaen and Shirakawa, [2] - existence; Yamada-Watanabe - uniqueness (see e.g., 
Rogers and Williams, p. 265 [10] or [6] p. 17 and Theorem 13.1) 
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possible on logarithmic scale. To this end we apply the Large Deviation 
Theory for asymptotic analysis of two families: 



} and I-t^Mt} 



where ^ 

Mt= [ -7)e-(^-^)^"dS,. (1.3) 
Jo 

For the second family, Large Deviation Principle (LDP) is well known. For 
the fist family, Freidlin-Wentzell's LDP, [5], is anticipated even though the 
diffusion parameter is only Holder continuous and singular at zero. For 
7 = ^, LDP is known from Donati-Martin et al.,[3], with the speed rate 
and the rate function of Freidlin-Wentzell's type with a corresponding 
modification: Jt{u) = ^ ( "'^/^"' )^'^{"t>0}'^^' ^ show that for 7 G 
(^,1) LDP is also valid with the speed rate and the rate function depending 
on 7. Combining both LDP's we obtain the following asymptotic result: 
there is a smooth nonnegative function , with Uq = 1 and absorbed at the 
time T, it^ = 0, such that for any smooth nonnegative function ut, with 
Uq = 1 and absorbed on the interval [0,T], ut = 0, 

lim lim — -r- — r- log P ( sup ~ < 5 



te[o,T] 



> lim lim — -— — r log P ( sup \xf^ — Uf\ < 5 

The latter inequality give us a motivation to consider uj' as the most likely 
path to absorbtion of the normed process xf. 

Note that calculations for P{tq < T) on logarithmic scale requires a non- 
standard technique. The set {tq < T} = {(xf )ig[o,T] £ D}, where 

D = {n G C[o,T] : = l;?^t = 'W6l(«)At, 6»{n)=mf{t:nt=0}<r}- 

D is closed in the uniform metric {g) in the space Cjqt] of continuous 
functions on [0, T]. Hence, the upper limit lim -7^57^^ log P (tq < T) is 

done according to the LDP technique. However, D has an empty inte- 
rior. This fact prevents us to use the LDP technique for the lower bound 
lim ^2(1-7) log P('^o ^ ^)- Nevertheless, we obtain this lower by using an 

inclusion {tq < T} ^ | < where Mx is defined in (1.3). The 

probability of the latter is easily computable, and gives a surprising result 

lim — i — - log \ -^—Mt < -1} = Ih^^ — rlogP(To<r). 

Of course, 

lim — ^ log f —^Mt < -i| < lim — . log P(ro < T) 



^For 7=5 see also [7] and Rouault [11] 
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which together the above establish the desired limit. This trick is of inde- 
pendent interest and might be useful for establishing LDP in other problems. 

2. Asymptotic of P(to < T) as K — > oc on logarithmic scale 

The random process Mt (see (1.3)) is a Gaussian martingale with the 
variation process {M)t = ^Mf: 

{M)t= r(j2(l-7)2e-2a-T)/^^ds. (2.1) 
Jo 



Theorem 2.1. For any T > 0, 

i^^^i=^^°^^("°^^) = -2(i^- 

Proof. To apply the Ito formula in a vicinity of tq, let us define a stopping 
time Tg = inf{t < T : xf" = e}, e > 0. Now, by Ito's formula, applied to 
(xfy-'f, t<TsAT,we find that 



)i-^ = 1 + Al - 7Mx^)'-^ds + [\l - 7)a 
Jo Jo 



dBs 



1 /■* ^2 1 

— / 7(1-7)^^71 — T^^TTT — ds (2.2) 

and, in turn, 

(r^ ^^^"^ -(l-7Mr.AT) 





In view of linig^o -^r^AT = -^roAT a.s. and the monotone convergence theo- 
rem 

0-2 g{l-7)MS 

in both sides of the above equality lim^^+o is applicable, that is, we have 

< (x^^^)i-^e-(i-^)'^(^°^^) 

f fj^ e{i-7)Ms 1 
+ / — ^7(1-7) -^ — ds =1 + —- M^AT- (2.3) 

(2.3) implies l + ^^MroAT > 0. If lj G {tq > T}, then 1 + ^^Mt{oj) > 0. 
In other words, {tq > T} C {1 + j^Mt > 0}, and so we obtain inclusion 



{to < T} D {j^Mt + K 0}. (2.4) 

It is well known that the families { j^i—y Mt } j^^^^ obeys LDP in the metric 
space (M, p) {p is the Euclidian metric) with the rate speed ^2(1-7) ^-i^d the 
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rate function I{v) = jITfY^- accordance with the large deviation theory, 



lim — i — ^logPfr < To) > - inf I(v) = . (2.5) 

^2(1-7) 6 V - u; _ ^.^^^^Q \ I 2{M)t 

A verification of the upper bound 
is more involved. We select a set 

D = {u e C[o,T] -^0 = l;Ut = e{u)=mf{t:ut=0}<T} 

which is closed in the uniform metric (g) related to the space Cjq^t] of con- 
tinuous functions on [0,T]. Obviously, {tq < T} C { (xf )ig[o,T] ^ D}, which 
suggests to find 

iSo ^ {i-f)teio,Tl € d) . (2.7) 

The most convenient tool to this asymptotic analysis is LDP for family 
{{^t^) t(z[Q T]} K^oo having the speed rate ^2(1-7) (0 and the rate function 
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Ut-llUtY dut=utdt, 

ut, r r"t-M"t i2„ 



00, otherwise 
(Theorem A.l). In accordance to the large deviation theory 

so that, it remains to prove 

inf Jt{u) 



ueD ' ' 2{M)t 
A minimization procedure of J{u) in u G D exclude from consideration 
functions ut with dut it. dt and Jq^'^^^'^ [ "^'t'"* ] = 00. This minimization 
is realized with a help of a specific deterministic control problem with a 
control action loj and a controlled process , being the solution of differential 
equation 

Ut = fjut + aujwt, t < e{u) A T (2.8) 

subject to the initial condition no = 1. Obviously, the function ut belongs 
to D. The pair (loj', 6{u*)), with u*{t) related to w*, is said to be optimal if 

/■e{n*)AT ne{u)AT 

/ iwtfdt < / w^dt 

Jo Jo 

for any pair {wt,ut) with J^^^^^"^ w^dt < 00. Technically, it is convenient to 

use the following change of variables: vt = u] ^ enables us to reduce the 
problem to a linear differential equation 

vt = lJi{l - l)vt + a{l - -i)wt] (2.9) 



instead of nonlinear one (2.8), subject to the initial condition vq = 1. We 

shall exploit also the property vaf^) , 

^ ^ 0, > 0. 

of equation (2.9), under the assumption 9{u) < T, implies: 



0, uq(u) =0 

^ ' The explicit solution 



= 7;e(„)e-''(i-^)* = [l + a(l-7) 



or, equivalently, the equality: 

1 



[0,e(n)AT) 



. , - -.-^^ (2.10) 

that, due to the Cauchy-Schwarz inequality, can be transformed into the 
inequality: 

2/i 2fi 



wfdt > 



> 



The choice of is conditioned by two requirements: 
1) (2.10) remains valid for wt replaced by 



Both requirements are satisfied for = —■^ -^_^_2ili—f)T e~^^"^~'^^*. Hence, 



Remark 1. ut 



(Mil 



Remark 2. The fact that the random variable Mt is gaussian with parame- 
ters (0, {M)t) and (2.4) yield for any K > 0,P{to <T) > P{Mt < -K^~'^). 

3. Most likely path to ruin of the normed process 

Since u\ = (f^)^/'^, where vl solves the differential equation 
v*t = /^(l - l)v*t + o-(l - l)w*t 



with -Wq = 1, we find that 



1 



, _ 2ii{M)t 

^ <72(l-7) 



<72(l-7) 



1/(1-7) 



1 



-2(l-7)Mi ■ 
-2{1-7)mT 



1/(1-7) 



On the other hand, in accordance with Theorem A.l for u* , we have 
1 



lim lim 



logP( sup \xf — < 6 

.t6[0,T] 



At the same time for any ti € D, Theorem A.l provides 
^— log PI 



-Jt{u*) 



hm lim , 
S^OK^cxi iv:2(l-7) 



logP( sup \xf — Ut\ < S 



-Jt{u) < -mf{u e D° < -Jt{u* 



Consequently, the function u* can be considered as the most likely path 
to ruin of the normed process on time interval [0, T]. 
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Appendix A. LDP for the family {(a;^)tg[o,T]}j<r^oo 

The family )t6[o,T]}x^oo is in Freidlin-Wentzell's framework [5]. In 
our setting, we take into account that the random process is absorbed 
at the stopping time tq, so that its paths belong to a subspace C^(J'^j(M+) 
of C[o^T](I^+) the space of continuous nonnegative functions ut = ^^^(u)) 
where 9{u) = inf{i <T:ut = 0}. The subspace C|j^^j(M4.) is closed in the 
uniform metric g and, consequently, it suffices to analyze the LDP in the 
metric space (C|'o^^](R+), g). The use of (C|'o^^](M+), g) instead of C[o,t](]R+) 
enables us to apply standard approach to LDP proof adding a few simplest 
details only. 

Theorem A.l. The family {{x^)t>o}K-^oo obeys LDP in the metric space 
(CjjJ'I^j (M+), £)) with the speed rate j^2(\~-y) ^'^^ '"'^^^ function 



1 / ( * ~ t\ /I+ dut=utdt, 

t [0,e{u)AT] 

oo, otherwise. 



Proof. The family { (xf )j>o}a'^oo is exponentially tight (see, e.g., theorems 
1.3 and 3.1, Liptser and Puhalskii, [8]), that is, 

lim lim^ ^^9n_^^ logPf sup > c) = -oo, (A.l) 



C^oo A-^0 i^2{l-7) ^ ^^jp 

l™n })^. P( sup Ix^+t - I > 7?) = -OO, (A.2) 

where r] is arbitrary number and is stopping time relative to a correspond- 
ing filtration. In (A.l), without loss generality xf might be replaced by 
{x^Y~"' what makes possible, in accordance with (2.2), to use the inequal- 
ity [xff'^i < e^^-^l^^'^[l + e(-i-''')^^(l - 7)0" ;|^], making the proof 
transparent. Due to (A.l), the condition from (A.2) can be replaced by an 
easy provable condition (here 2lc = { supj<j' xf < C}): 

lim lim sup — rr- — r log pf sup \x^,f — x^\ > ri, = —00, V C > 0. 
For e(u) > T, the proof of local LDP: 

does not different from standard one and is omitted. The case of Jt{u) = 
— 00, including 7^ 1, dut ^ dt, is analyzed in a standard way and is 
omitted too. 

The analysis of "uq = 1, dut = iitdt, /^^^"^^^ < oo> ^(^) < ^" 

is based on the following result. 

Proposition A.l. [Dupuis, Ellis [4], A. 6. 3.] For any absolutely continuous 
function u = ('Ui)tg[o,T]j mapping [0,T] into and any a € M 

I (ut=a\ dt = 0. 
Iwt^O/ 



Local LDP upper bound. Set = ^ V ut and notice that 0{u"') > 
T. Moreover, Uq = 1, duf = iifdt and, due to Proposition A.l, it" = 

^t^{„,>i}ds and also /g^^" ( "^"n^^" )^^'^ ^ 

Since r" = inf{i : < -} ^ ^(^^)) ^ oo, we find that 



hm lim — -7- — r- log P ( sup I xf — uA < 6 

9{u)/\T 



n^oo 2 



Local LDP lower bound. With (/> > (5 > 0, write 
I sup |xf -<| < 5] 

= 1 sup |xf - n"| < (5[ n I sup |n"-nt|<0[ 
tG[o,T] te[o,T] ^ 

Ml sup |xf - u"| < (ji n I sup \vJl-ut\>4>\ 

c| sup - ut| < </> + (5[ n I sup |u"-nt|>(/>[ 

c| sup \xf - ut\ < 2(j)\ \J \ sup |n"-Uf|>(^[. 
tG[o,T] J ^ ie[o,T] 

For fixed cf), there exits a number > ^ such that for any n > the set 
I supjg[o,T] 1^? — I > i;^} = 0- Therefore, for sufficiently large numbers n, 

P( sup |xf - uti < 2(/>) > P( sup |xf - uj"! < 5). 
Hence and by Proposition A.l, a chain of lower bounds holds, 

> lim lim ^ , logPf sup \xf — u^l < 6 



1 /u,-^n,x2 



2(t2 
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providing 



IT^fTTvrlogPf sup \xf-ut\< 



te[o,T] 



□ 
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